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Abstract 

We present an analytical result for the supercurrent across a 
superconductor/quantum-dot /superconductor junction. By converting the 
current integration into a special contour integral, we can express the cur- 
rent as a sum of the residues of poles. These poles are real and give a natural 
definition of the Andreev bound states. We also use the exact result to explain 
some features of the supercurrent transport behavior. 
PACS numbers: 74.50.+r, 73.40.Gk, 73.20.-b, 73.63.Kv. 

Introduction. The mesoscopic superconductor/quantum-dot /superconductor (S-QD-S) 
junction is a typical structure to study the phase coherent current transport in mesoscopic 
hybrid systems [|TJ . When using Keldysh Green function technique [0 to evaluate the current 
across the S-QD-S junction, the usual method demands a numerical integration and faces 
the problem of explaining some features of the current transport behavior ||. In this paper, 
we employ a new method to analytically compute the current across the S-QD-S junction in 
the absence of voltage bias, and use the analytical results to explain some features of current 
transport properties. 

Model and Hamiltonian. In the system under consideration, a quantum dot defined in a 
2-dimensional electron gas (2DEG) is coupled to two BCS superconducting leads. We model 
this system by the Hamiltonian 

H = H L + H d + H R + H T , (1) 
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are the isolated (unperturbed) Hamiltonians of the left superconductor, the quantum dot, 
and the right superconductor, respectively. 

We consider only a quantum dot with a negligible intra-dot Coulomb interaction. The 
reason is that we want to focus on the superconducting proximity effect on the QD, which 
is due to Andreev reflections; whereas a large charging effect caused by intra-dot Coulomb 
repulsion suppresses the transport mediated by Andreev process . Study on the competi- 
tion of a medium intra-dot Coulomb interaction with Andreev reflections in S-QD-S junction 
will be carried on in the coming work. Furthermore, we consider only one single energy level 
inside the QD. The generalization to QDs with several energy levels is straightforward. 

The tunneling term 

Ht = Y\ L Aj° + c.c] + Y,l R P b l*da + c.c.} , (5) 
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describes the electron-transfer between the QD and the leads. 

Current Formula. The current flowing across the junction is || 
2er r „ r+™ du 
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where g^, G are the Keldysh Green functions of the left superconductor and the QD in 
Nambu space M. The retarded Green function of the superconductor is 
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solving the matrix Dyson equation, we obtain the retarded Green function of the QD: 



where p L / R {uo) = , m is the density of states (DOS) of a BCS superconductor. After 
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is the unperturbed retarded Green function of the QD. The QD's self-energy 
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describes the proximity effect of superconductors on the QD; here we have used the wide 
band-width approximation pj. D(uj) in Eq. [8] is the determinant of the matrix (g r_1 — S r ): 



D(uj) = [g r (a;)" 1 - 57 (u,) n ][<f (a;)" 1 - X7 (u;) 22 ] - E r 57 ( w ) 
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Choice of Integration Contour. The supercurrent formula Eq. [6] can be expressed as 

he = ~ Re{ f°° ^[G d (u) 12 g< (u) 21 + G< (uj) 12 g a L (u) 21 ]} . (12) 

h J-00 Z7T 

In order to avoid calculating this integration numerically, we shall transform it into a contour 
integral. First, using the fluctuation-dissipation theorem G < (uj) = [G a (uj) — G r (u)]f(uj), 
where G stands for any Green function in the equilibrium state and f(uj) is the Fermi 
distribution function, we have 
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Then we make a change of the integration variable u —>■ uj — i0 + for G r d (ui) l2 g r L 
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and uj — > uj + i0 + for G a d (uS) X2 g a L (uj) 21 f(u>) to divide I sc into two integrals along different 
paths: 

ApV /-+oo+i0+ rjrj . . . 
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Since the above two integrands are actually the same, we define a new integrand equal to 
them: 

J{u) = G d (u - i 0+) i2 g r L (u - z0 + ) 2i f(u - t0 + ) 

= G a d (u + *0+) i2 gl (u + *0+) 2i /( W + *0+) (15) 

Thus the current integration can be rewritten into a contour integral in the complex-u; plane: 

4er^ f +oo+io+ ,+oc-io+ d u 4 e r r dcu 

Isc = - r MU -/ W(w =-=-Ite/ja;— . (16) 

The integration path C = f^t^+ + J7^Z^q+ is a close contour which lies infinitely close to 
the real u axis. The integral can then be evaluated analytically using Cauchy's theorem. 

This choice of the integration contour has the advantage of leaving all (infinitely many) 
the poles at uj — ^ 2n + 1 ) 7r arising from f(u) outside the contour, thus avoid the calculation 
of the summation over all the Matsubara frequencies (u = i ^ 2n+ ^ ). 

Poles and definition of Andreev bound state. Since the current is proportional to the 
sum of the residues of the integrand J{oS) inside the contour C, only real poles of J{oS) 
contribute. For J(u) = G r d (cu — i0 + ) l2 (?£ (uj — ^0 + ) 21 f(co), there are only two pairs of real 
poles: uj = ±A and uj = ±e*(|e*| < A). The poles at u — ±A originate from the singularities 
of the DOS of the BCS superconductors. 

The poles at ui — ±e* come from the QD's Green function and are two real roots of the 
equation D(uo) = inside the gap ( |e*| < A ). These two poles give the positions of the 
quasi-particle states of the QD inside the gap. The fact that they are real indicates that 
these two quasi-particle states are exactly bound states, usually called the Andreev bound 
states (ABS). 

Before coupling to the superconductors, the unperturbed QD has a bound state at e 
for the electron spectrum and one at — e for the hole spectrum. After coupled to the 
superconductors, through the electron-transfer with the two superconductors via Andreev 
process, the electron part and the hole part of the QD's Green function are coupled together, 
and give rise to significant modifications to the QD's energy spectrum: the original bound 
state at e for electron (or — e for hole) renormalizes into two symmetrical bound states at 
±e* for both electron and hole spectrum. 

To express e* as the algebraic function of e , </> and T, we transform the equation D(uj) — 
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into a quartic equation of x = e* 2 : 



x 4 - (2A 2 - 4r 2 )x 3 + (A A + 2A 2 5 2 - 4r 2 A 2 )x 2 + (2A 2 A 2 B 2 )x + (A 2 B 2 ) 2 = 0; 

A 2 = A 2 + e 2 + T 2 , B 2 = e 2 + T 2 cos 2 \ ; x = e* 2 < A 2 . (17) 

This quartic equation of x has only one positive real solution xq less than A 2 : ±e* = ±y/xQ. 
The quartic equation is easy to solve algebraically; but we omit the lengthy solution. The 
dependence of e* on e , (f> and Y is shown graphically in ||. 

Residues and analytical result of the current. The current is obtained by adding up the 
residues of the four poles: 

AeT 4 

I sc = — ImEResjy], (18) 
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where o>j = ±e*, ±A. The residues of J at e* and — e* are : 

Res. =±e , J(u>) = Jim [ "~ (±e V (w - ^0 + ) i2 ^ (w - ^0 + ) 2i /(«)], (19) 
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The residues at ±A turn out to be zero since 

D(±A) = oo, -^1^ = 0. (20) 

Thus the supercurrent is simply 

e r 2 A 2 sin^tanh^f 
Lc = ~h~ A 2 - e* 2 e *[(A 2 - e* 2 ) + 7 ^ ? f(2A 2 - e* 2 ) + ^^A 2 sin 2 f] ' (2 " 

where = — 0l is the phase difference of the two superconductors. 

Discussion. In addition to the splitting of the original energy level eo into the two 
ABS, another modification to the QD's energy spectrum is that the QD begins to have 
a continuum energy spectrum outside the gap. In fact, one way to calculate the current 
integral ( I sc = -if Re {J^^[G r d (u;) 12 gf(u;) 21 + G< (a,) 12 g a L (u,) 21 ]} ) is to divide the 
QD's spectrum into the continuum part ( \u\ > A ) and the discrete part ( \u\ < A ), and 
then evaluate their contributions to the current separately. However, this method involves 
a numerical integration and faces the problem of explaining why the peaks at u — ±A in 



5 



Icontinuum vs. €q curve and in Idiscrete vs. 6q curve exactly cancel out when the two parts add 
up to give the total current, see Fig.l. 

In contrast, from our approach, this exact cancellation can be shown naturally. The cur- 
rent is obtained by adding the residues of the four poles. Since the residues at ± A vanish, the 
total current curve has no peaks arising from the DOS singularities of the superconductors. 

Now only residues at ±e* contribute to the total current and give rise to a sharp peak in 
/ vs. eo curve at eo = 0. The reason for the appearance of the peak at eo = is that when 
eo = 0, the upper and the lower ABS have same phase. As eo departs from the Fermi level, 
the phase of the upper and lower ABS begins to differ and this difference severely reduces 
the supercurrent across the junction. Fig. 2 shows the central peaks at different T, the peaks 
broaden when T increases. 

A direct application of this result is a way to align the energy level of the QD to the 
Fermi level of the superconductor in the actual experiment: when changing the gate voltage, 
the state with the maximum supercurrent is at e = 0. 

Conclusion. In summary, we have calculated analytically the supercurrent across an S- 
QD-S junction. The central results are Eq.[18] and Eq.[21]. They give the analytical result 
of both the positions of ABS and the supercurrent. We show two additional advantages 
of our method: (1) It presents a natural definition of the Andreev bound states inside the 
QD. (2) It explains the disappearance of the peaks at the gap edge in / vs. e curve. The 

method devised in the paper can also be used in computing other equilibrium properties of 
the superconductor-normal hybrid systems. 
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Science Foundation of China under the Grant No. 10074001. 
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FIGURE CAPTIONS 



Fig. 1 Supercurrent vs. eo when T = 0.1 and (f) — | . The dash line is from discrete spectrum, 
and the dot line from continuum spectrum. The solid line is total current. Note that 
the continuum spectrum contributes negative current. 

Fig. 2 Supercurrent vs. e with <p = | for different T. Different curves are for T = 0.1, 0.3, 
1.0, respectively. Note the broadening effect as T increases. 
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Fig.1. 
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